Elliptic flow of (p meson a sensitive probe for onset of de-confinement transition 

in high energy heavy-ion collisions 
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Elliptic flow (vi) of <f> meson is shown to be a sensitive probe of the partonic collectivity using 
A Multi Phase Transport (AMPT) model. Within the ambit of the AMPT model with partonic 
interactions (string melting version), the <j> meson vi at midrapidity is found to have negligible 
contribution from hadronic interactions for Au + Au collisions at ^/s KN = 200 GeV. Changing the 
hadron cascade time in the model calculations does not change the cj> meson V2, while it reduces 
the i>2 of proton at low transverse momentum (pr)- These observations indicate that a substantial 
reduction in <j> mesons V2 as a function of colliding beam energy would suggest the dominance of 
hadronic interactions over partonic interactions. 

PACS numbers: 25.75.Ld 



I. INTRODUCTION 

One of the current focus of the high energy heavy- 
ion collision experiments is to study the various 
aspects of the QCD phase diagram [l|, After 
observing the clear signatures of the formation of 
strongly interacting quark-gluon plasma (QGP) at 
the top RHIC energies attempts are being 

made to vary the colliding beam energy and search 
for the transition region between the partonic and/or 
hadronic dominant interactions in the phase dia- 
gram. This is one of the goals of the RHIC Beam 
Energy Scan (BES) program 

In Au + Au collisions at y/s^" = 200 GeV the <f> 
meson production has played a crucial role to estab- 
lish the formation of the partonic matter Q . Num- 
ber of constituent quark scaling of the elliptic flow 
of (f> meson [Toj , enhancement in yield of <f> meson in 
Au + Au collisions relative to p + p collisions [ll[ 
and the ratio of yield of the baryon to the yield of 
4> meson as a function of pt fl0| have been the key 
measurements. In general, multi-strange hadrons 
are understood [l2[ to have a smaller hadronic in- 
teraction cross sections. As a result, their momen- 
tum distributions, and particularly the momentum 
dependence of the azimuthal anisotropy [l(|, are 
expected to be primarily controlled by the early 
partonic interactions in high-energy nuclear colli- 
sions. In order to demonstrate that the multi- 
strange hadrons, are a clean tool for phase identifica- 
tion, we have performed model simulations to study 
the (f> meson production. The paper focuses on us- 
ing the AMPT model for the <j) meson production, 
leading to the establishment of the <f> as a key ob- 
servable for studying the onset of the de-confinement 
transition. 



Further, with high statistics data being collected 
in high energy heavy-ion collisions, it is now possible 
to have high precision measurements of v% for various 
produced hadrons. One possibility that opens up 
is to study the effect of the late stage hadronic re- 
scattering on vi at low px- Initial simulations using 
a hybrid model showed that the usual mass ordering 
trend of V2(4>) < ^2(p) will get reversed due to the 
late stage hadronic re-scattering [l3|]. In this paper 
we also study this aspect using the AMPT model. 

The paper is organized as follows. In the 
section HH we will briefly introduce the AMPT 
model [Tj] used in this study. In the section [TTTT the 
results from the AMPT model calculations regard- 
ing V2 of <f> mesons and protons for various config- 
urations (different parton-parton cross section and 
hadronic cascade time) are presented. Finally in 
section IIV1 we summarize our findings and present 
a short discussion on the implications of this work 
to the current experimental measurements in high 
energy heavy-ion collisions. 



II. MODEL CALCULATIONS 

A. AMPT 

The AMPT QJ] model used for the calculations 
presented in this paper has four main stages: the 
initial conditions, partonic interactions, the conver- 
sion from the partonic to the hadronic matter, and 
hadronic interactions. The initial conditions are ob- 
tained from the HIJING model [r|. Scatterings 
among partons are modeled by Zhang's parton cas- 
cade (ZPC) QJ], it includes only two-body scatter- 
ings with cross sections obtained from the pQCD 
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with screening masses. Some of the results pre- 
sented are by varying the parton-parton scattering 
cross sections within 3 mb to 14 mb. The AMPT 
model with string melting [TtJ leads to hadron for- 
mation using a quark coalescence model. The sub- 
sequent hadronic matter interaction is described by 
a hadronic cascade, which is based on A Relativis- 
tic Transport (ART) model [ill- The termination 
time of the hadronic cascade is varied in this pa- 
per from 0.6 fm/c to 30 fm/c to study the effect of 
the hadronic re-scattering on the observables pre- 
sented. More detailed discussions regarding the 
AMPT model can be found in Ref. [lj]. I 11 this 
study, approximately one million events for each 
configuration (different cross section and hadronic 
cascade time) were generated for Au + Au 0-80% 
minimum bias collisions at ^/s NN — 200 GeV. All 
results presented are for the rapidity range ± 1.0 
unit. 



C. Elliptic flow 

The anisotropic elliptic flow parameter presented 
in this paper is defined as the 2 nd Fourier coefficient 
V2 of the particle distributions in emission azimuthal 
angle (<p) with respect to the reaction plane angle 
(^) [2(j, and can be written as 

dN 

— oc l + 2i J2 cos(2(0-tf)). (1) 
acp 

For a given rapidity window the second coefficient is 

2 2 

« 2 = <cob(2(0-*))) = (^^), (2) 

Px+Py 

where p x and p y are the x and y component of the 
particle momenta. In the AMPT model the ^ is 
along the x-axis. 



III. RESULTS 



B. 



mesons production 



The string melting version of the AMPT used in 
this paper, produces <j> meson using a quark coa- 
lescence model in the partonic stage pjj]. In the 
hadronic stage, the AMPT model includes the fol- 
lowing reactions associated with the <f> meson. In- 
elastic scatterings in baryon-baryon channels in- 
cludes (NAN*)(NAN*) -> <f>NN, those in the 
meson-baryon channel includes (irp)(NAN*) <H- 
(j)(NAN*) and K(AE) O <j>N. The <j> meson scat- 
terings with other hadrons included in the model 
are <j){-7tpLo) <-> (KK*)(KK*), and <j){KK*) O 
(tt pu))(K K*) . The cross section for the elastic scat- 
tering of the 4> meson with a nucleon is set to 8 mb 
while the <f> meson elastic cross section with a me- 
son is set to 5 mb. For other details and specifically 
those related to inelastic scattering cross section can 
be found in 

For comparison with another transport based cal- 
culation, we also present V2 for <f> meson from 
UrQMD (Ultra relativistic Quantum Molecular Dy- 
namics) model [lj| for Au+Au collisions at ^/snn 
= 200 GeV. It is based on a microscopic transport 
theory where the phase space description of the re- 
actions are important. It allows for the propagation 
of all hadrons on classical trajectories in combina- 
tion with stochastic binary scattering, color string 
formation and resonance decay. It incorporates 
baryon-baryon, meson-baryon and meson-meson in- 
teractions, the collisional term includes more than 
50 baryon species and 45 meson species. 
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FIG. 1: (Color online) <j> meson V2 for Au + Au mini- 
mum bias (0-80%) collisions at midrapidity (± 1.0) at 
y/s^ = 200 GeV from STAR experiment at RHIC H(} 
compared to the corresponding AMPT and UrQMD 
model calculations. The AMPT model calculations are 
shown for three different parton-parton interaction cross 
sections of 3, 10 and 14 mb. The errors shown are sta- 
tistical. 

Figure [T] shows the comparison of elliptic flow of 
4> meson in 0-80% minimum bias Au + Au collisions 
at midrapidity for y/s NN = 200 GeV measured by 
STAR experiment at RHIC [uj with the correspond- 
ing results from the AMPT model for three different 
values of parton-parton corss sections of 3 (magenta 
solid circle), 10 (black open square), and 14 mb (red 
open cross). Although it is generally believed that 
perturbative QCD cross section is about 3 mb (2lj . 
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in order to explain the measurements for pt > 1 
GeV/c a parton-parton cross section between 10-14 
mb is found to be required. Thus the generation 
of this substantial vi for cj> mesons as observed in 
the experiments requires a significantly large parton 
interaction cross section than obtained from pQCD 
calculations (~ 3 mb) 0, [H [H- Also shown in 
the figure is the corresponding vi results from the 
UrQMD model [H[. The UrQMD model results 
(which docs not include any partonic interactions) 
gives substantially smaller value of vi compared to 
the experimental data. 

Within the framework of hydrodynamics, the col- 
lectivity reflected through the V2 distributions are 
caused by the pressure gradients. Due to the effect 
of self quenching, one expects that the development 
of vi is dominantly from the early stages of the col- 
lision [22I, [23| . From the comparison shown in Fig. 
1, we conclude that collectivity in the multi-strange 
hadron meson has been developed in the early par- 
tonic interactions in high-energy nuclear collisions at 
RHIC 0|. We now proceed to investigate the con- 
tributions to </> meson vi in AMPT separately from 
partonic and hadronic interactions. 

A. (f> meson v-z from partonic interactions 



sus pt from AMPT model for parton-parton cross 
section of 10 mb (red solid circles) and results with- 
out any parton-parton interaction (blue solid square, 
obtained by setting the parton-parton cross section 
value to mb). The hadronic cascade time is 30 fm/c 
for both the cases. The <j> meson V2 is consistent with 
zero in absence of parton-parton interactions. The 
panel (c) shows the difference between the two re- 
sults, indicating that almost all the <f> meson V2 is 
generated via the partonic interactions. 

Figure [2] (b) and (d) reinforces these observations 
by presenting the <p meson V2 for the same system be- 
fore (solid blue squares) and after (solid red circles) 
the relativistic transport calculations for hadrons. 
The results are similar between the two cases as seen 
from the px dependence of <fr meson V2 in panel (b), 
the difference between the two cases shown in panel 
(d) is consistent with zero. These results indicate 
that partonic interactions are essential for generat- 
ing V2 of c/> mesons and contributions from hadronic 
interactions are minimal, within the context of the 
AMPT calculation. 



B. Effect of hadronic re-scattering 
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FIG. 2: (Color online) (j> meson vi for Au + Au mini- 
mum bias (0-80%) collisions at midrapidity (± 1.0) at 
y/s^ = 200 GeV from the AMPT model. Panels (a) 
and (b) shows the results as a function of pr for parton- 
parton interaction cross section of and 10 mb and cal- 
culations before and after relativistic transport calcula- 
tions for hadrons, respectively. The lower panels (c) and 
(d) shows the difference in V2 shown in panels (a) and 
(b), respectively. The errors shown are statistical. 

Figure [2] (a) shows the <f> meson V2 for mini- 
mum bias Au + Au collisions at midrapidity ver- 
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FIG. 3: (Color online) (a) «2 of protons as a function of 
p T for Au + Au 0-80% collisions at y/sj^ = 200 GeV 
from AMPT model at midrapidity. The results are 
shown for a parton-parton cross section of 10 mb and 
three different values of hadronic cascade time periods, 
(b) The same plot as (a) for the <j> mesons, (c) Ratio of V2 
of protons for hadron cascade time of 0.6 fm/c to corre- 
sponding 1)2 for time periods of 15 and 30 fm/c, and (d) 
same as in (c) for the <f> mesons. The error bars shown 
are statistical. 

To further study the effect of hadronic interac- 
tions, the model simulations were carried out for 
the Au + Au minimum bias collisions with parton- 
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partem interaction cross section fixed to be 10 mb 
and varying the hadronic cascade time from 0.6 
fm/c to 30 fm/c. Higher value of hadronic cascade 
time reflects larger hadronic re-scatterings. We have 
checked that for RHIC energies going to even longer 
time duration does not contribute any further to the 
results presented. The calculations arc carried 
out for both <f> meson and proton. We chose protons 
mainly for two reasons: (a) as a hadron, it has a 
mass similar to that of the <fi meson and (b) con- 
trary to that of the <fi meson, it has larger hadronic 
interaction cross sections. 

Figure [3] (a) shows the V2 of protons versus pt 
for 0-80% Au + Au collisions at ^/s^ = 200 GeV 
from the AMPT model with parton-parton inter- 
action cross section of 10 mb and three different 
values of hadronic cascade time of 0.6 fm/c (red 
solid circle), 15 fm/c (black open cross) and 30 fm/c 
(blue solid square). With increase in hadron cas- 
cade time, which reflects increasing contributions 
from hadronic interactions, the proton vi decreases 
at lower px- Implying the development of the collec- 
tive expansion in the hadronic period. This is more 
clearly illustrated in the panel (c) of the figure, which 
shows the ratio of the proton V2 for the hadron cas- 
cade time of 0.6 fm/c to the corresponding v 2 values 
for time periods of 15 (open crosses) and 30 (solid 
squares) fm/c. Figure [3] (b) and (d) shows the cor- 
responding results for <fi mesons. In marked contrast 
to the case for protons, the <f> meson vi remains unaf- 
fected by the hadronic interactions, indicating that 
vi is solely generated due to the partonic interac- 
tions in these model calculations. 

Figure 2] shows the effect of hadron re-scattering 
on the mean values of the transverse momentum dis- 
tributions for proton and <j) meson. The average 
transverse momentum < pt > for protons increases 
as the hadronic cascade time increases, where as for 
cj) mesons this change is much smaller. The results 
are shown for two different values of parton-parton 
cross section values of 3 and 10 mb. For both cases, 
the trends are identical, suggesting that the increase 
in value of < pt > with hadron cascade time is 
dominantly due to hadronic interactions. The mag- 
nitude of < pt > though is decided by both partonic 
and hadronic interactions taken together. Hence the 
< Pt > for calculations with parton-parton cross 
section of 3 mb is lower compared to those with 
parton-parton cross section of 10 mb. 



IV. SUMMARY 

In summary, we have studied the generation of 
elliptic flow in AMPT model for <f> mesons at midra- 
pidity for 0-80% centrality Au 4- Au collisions at 
y/s NN — 200 GeV. The large value of V2 for meson 
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FIG. 4: (Color online) Average transverse momentum 
(< pr >) for proton and cj> meson produced in 0-80% 
Au + Au collisions at midrapidity for Js NN = 200 GeV 
from AMPT model. The results are shown as a function 
of hadron cascade time and for two values of parton- 
parton cross sections. 



from the model which is comparable to the corre- 
sponding measurements at RHIC is due to the par- 
tonic interactions. Hadronic interactions as mod- 
elled in AMPT arc unable to generate any vi for 
4> meson. Increasing the hadronic re-scattering by 
increasing the hadron cascade time period does not 
change the conclusions. In contrast to (f> meson, the 
i>2 of protons decreases with increase in hadronic re- 
scattering effects at low transverse momentum. This 
effect is also seen in the average transverse momen- 
tum of protons which increases with increase in the 
time duration of the hadron cascade. 

We have identified the multi-strange hadron </> me- 
son as a clean tool for studying partonic interactions 
in high-energy nuclear collisions. Any significant 
change in <p meson collectivity will signal the pos- 
sible change of in the nature of the medium formed 
in high energy heavy-ion collisions. Zero value of the 
i>2 for the 4> meson would clearly indicate that the 
system formed in heavy-ion collisions did not make 
the de-confinement transition. Hence </> meson vi is 
an ideal probe for the search of the phase boundary 
in the QCD phase diagram. 

In addition our study shows that proton vi at low 
Pt decreases with increase in hadron cascade time. 
This rcflcctes the effect of late stage hadronic re- 
scattering. Whereas for the <fi meson the vi remains 
unaffected by late stage hadronic re-scattering pro- 
cess. In experiments it has been observed that at low 
Pt a distinct mass ordering of vi is followed. Heavier 
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particles have smaller [24J . If re-scattering is suffi- 
ciently large the proton vi at low pt could be smaller 
than the corresponding V2 for the <j> meson. The key 
fact being that the V2(4>) remains unchanged with 
increase in hadronic re-scattering, while the K2O?) 
decreases at low px- This can then be considered to 
be a signature of hadronic re-scattering and can be 



probed using high event statistics data sets in the 
high energy heavy-ion collision experiments. 
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